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Book Shelf - TIME 2008 


D-N-L#69 


[1] Rainer Wonisch , Parabeln und Co. erforschen mit TI-Nspire CAS (Software) 
bk teachware SR-58, ISBN 978-3-901769-79-5 

[2] Rainer Wonisch , Parabeln und Co. erforschen mit TI-Nspire CAS (Handheld) 
bk teachware SR-59, ISBN 978-3-901769-80-1 

{Inhalt von [1] und [2]: Geraden im Koordinatensystem, Parabeln im Koordinatensystem, Null- 
stellen von Parabeln, Anwendungen ) 

[3] Josef Bohm , Programmieren mit dem TI-Nspire CAS 
bk teachware SR-60, ISBN 978-3-901769-81-8 

{Inhalt: Anstelle einer Einfuhrung: die Finanzen, Uber eine Sutra zur ersten Bibliothek, Das 
Chaos-Spiel, Eine Bibliothek fur Dreiecke, Wir haben ein CAS: ein Programm fur Extrem- 
wertaufgaben, Trainingsprogramme fur Grundfertigkeiten, Wir iiberwinden die Regression 4. 
Grades, Mit TI-Nspire zum Lottohaupttreffer) 

[4] Pat und Carl Leinbachl Josef Bohm , Forensische Mathematik fur den Unterricht 
bk teachware SR-57, ISBN 978-3-901769-78-8 

{Inhalt: Die Formel von Sherlock Holmes, Kriminologie mit dem Wetterhauschen, Eine Buch- 
stabensuppe als letzte Nachricht, Todliche Geschwindigkeit, Die Botschaft der Bremsspuren, 
Blutspritzer an der Wand, Mit der Matrix zu den Tatverdachtigen, Knochenreste vom Bruder- 
krieg, Stumme Zeugen: Finger abdrucke, Ein GPS - nicht nur fur das Auto, Die Aufgaben eines 
Coroners) 

If you know about CAS-related publications then please inform me. Books dealing 
with the use of technology in math education in general are also interesting. Josef 


TIME 2008 Call for Papers 

22. September- 26. September 2008 

http : / / time . tut . ac . za 


• 30 April 2008: Submission of abstracts on or before this date 

• 31 May 2008: Acceptance of abstract notification by this date 

Abstracts can be submitted as e-mail attachments to Josef Bohm ( nojo.boehm@pgv.at ) or 
Michel Beaudin ( Michel.Beaudin@etsmtl.ca ) for the Derive/Nspire part of the conference 
and for the ACDCA part to Bernhard Kutzler ( b@kutzler.com ) or Vlasta Kokol-Voljc 
(vlasta.kokol@uni-mb.si ). 

The deadline for submission is 30 April 2008 with acceptance of notification by 31 May 
2008. 


Download all DNL-Derive- and Tl-files from 

http : / / www . austromath . at/ dug/ 
http : // www . bk- teachware . com/ main . asp?session=375059 
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Liebe DUG-Mitglieder, 

Mit dieser 1. Ausgabe des DNL im Jahr 2008 
mochte ich Si e nochmals recht herzlich zur Teil- 
nahme an der TIME 2008 einladen. Die Einreich- 
frist fiir einen Beitrag endet mit 30. April 2008. 
Vielleicht konnen Sie diese Konferenz auch mit 
einem weiteren Aufenthalt in Siidafrika verbinden. 

Nun zu diesem DNL#69: er ist wieder eine Mi- 
schung aus unterschiedlichen Einsatzgebieten wie 
Statistik (ANOVA), 3D-Geometrie (Tetraeder), 
Schulmathematik und Einstieg in die Hochschulma- 
thematik (Devoir #1). Diesen letzten Artikel (von 
Michel Beaudin) habe ich gerne in seiner Original- 
fassung belassen, zumal auch die Losungen beige- 
fugt sind. 

Ich wollte eigentlich auch noch einen Beitrag zu 
den FuBpunktkurven und -flachen liefern, aber - 
wie so oft - muss ich diesen fur den nachsten DNL 
aufheben. Als kleinen Vorgeschmack konnen Sie 
sich auf den Seiten 14 und 36 holen. 

Wenn Sie sich auch fur schone und interessante 
3D-Grafiken interessieren, habe ich einen Hinweis 
far Sie: das ausgezeichnete Programm Surfer kann 
frei bezogen werden. Sehen Sie bitte die entspre- 
chende URL im User Forum. (Einige damit erzeug- 
te Grafiken finden Sie in diesem DNL.) 

Von unserer Homepage konnen Sie neben den 
DNLs den aktualisierten Index fur DNL#1 - #16 
und nun neu auch den Index fiir DNL#53 - #68 
finden. Daneben habe ich von alien downloadbaren 
DNLs die Titelseiten zusammengestellt. Damit 
kann auch so nach Inhalten geblattert werden. 

In den letzten Marztagen gab es wieder ein Treffen 
mit Damen und Herren aus dem TI-Nspire- 
Entwicklungsteam. Es werden groBe Anstrengun- 
gen untemommen, diese Software auch den Wiin- 
schen und Erwartungen der DERIVE- und Voyage 
200 (TI-89/92) User anzupassen. Es wurde auch die 
Gelegenheit wahrgenommen, insbesondere Anre- 
gungen, Forderungen und Empfehlungen hinsicht- 
lich eines sinnvollen didaktischen Einsatzes eines 
CAS weiterzugeben. Viele Probleme in Bezug auf 
Beniitzerfreundlichkeit konnten offen diskutiert 
werden. Die DUG ist nach wie vor ein wichtiger 
Gesprachspartner fiir TI. 

Bei dieser Gelegenheit darf ich die Nspire- 
Anwender (Software und Handheld) wieder um 
Beitrage, Hinweise, Tipps und Tricks bitten. Natiir- 
lich sind auch weiterin alle anderen Beitrage sehr 
willkommen. Auf der Info-Seite konnen Sie einige 
Neuerscheinungen zum TI-Nspire finden. Auch auf 
der ACDCA-Homepage ( www.acdca.ac.at) gibt es 
Beitrage zum Programmieren mit dem TI-Nspire. 

Mit einem nochmaligen Hinweis auf die TIME 
2008 verbleibe ich Euer 


Dear Dug-Members, 

At the beginning of this 1 st issue 2008 I’d like to 
invite you once more participating at TIME 2008. 
Deadline for submitting an abstract is 30 April. 
Maybe that you can combine this conference with a 
holiday in South Africa. 

Some comments on this DNL#69: it is again a 
mixture of various applications as statistics 
(ANOVA), 3D-Geometry (Tetrahedron), secondary 
school math (From a Spreadsheet to a CAS) and 
first steps in higher mathematics (Devoir #1). I left 
the last article in its original version (French), the 
more so as it contains the solutions. 

I intended to deliver a contribution on pedalcurves 
and -surfaces but - as so often - I had to postpone 
it for the next DNL because of lack of space. On 
Pages 16 and 36 are shown mouthwatering exam- 
ples . 

If you are interested in beautiful and interesting 
3D-graphics then I have a tipp for you: download 
for free the great program Surfer. You can find the 
respective URL(s) in the User Forum. (Some Surfer 
produced graphs are presented in this DNL.) 

You are invited to download from our homepage 
the updated index for DNLs#l - #16 and new also 
the index for DNLs#53 - #68. Besides this I col- 
lected all frontpages of the said DNLs in two pa- 
pers. So you can easily browse the contents of all 
downloadable DNLs. 

In the last days of March we had a 3 -days meeting 
with people of the Tl-development team. There are 
big efforts to accommodate the CAS software to the 
wishes and expectations of the DERIVE- and the 
TI-89/92/V200 users. We took also the opportunity 
to pass a bundle of ideas, demands and recommen- 
dations with respect to a meaningful didactical use 
of a CAS. Many user friendliness problems were 
discussed very open. The DUG proved again to be 
an important interlocutor for TI. 

At this occasion Ed like to ask the Nspire-users 
(software and handheld) to submit contributions, 
comments, tipps and tricks. Of course, all other 
contributions are very welcome as ever. On our 
info page you can find some new publications for 
the TI-Nspire. You can also download contribu- 
tions (in German) from the ACDCA-homepage 
( www.acdca.ac . at ) . 

Repeating my reminder on TIME 2008 I remain 
your 
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The DERIVE-NEWSLETTER is the Bulletin 
of the DERIVE & CAS- 77 User Group. It is 
published at least four times a year with a 
contents of 40 pages minimum. The goals 
of the DNL are to enable the exchange of 
experiences made with DERIVE, 77-CAS 
and other CAS as well to create a group to 
discuss the possibilities of new methodical 
and didactical manners in teaching mathe- 
matics. 


Editor: Mag. Josef Bohm 
D'Lust 1, A-3042 Wurmla 
Austria 

Phone/FAX: 43-(0)2275/8207 
e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. 
It must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS-77 Newsletter will be. 

Next issue: June 2008 

Deadline 15 May 2008 


Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 

Wonderful World of Pedal Curves, J. Bohm 

Tools for 3D-Problems, P. Luke-Rosendahl, GER 

Financial Mathematics 4, M. R. Phillips 

Hill-Encription, J. Bohm 

Farey Sequences on the Tl, M. Lesmes-Acosta, COL 
Simulating a Graphing Calculator in DERIVE, J. Bohm 
Henon & Co, J. Bohm 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Liike-Rosendahl, GER 

Overcoming Branch & Bound by Simulation, J. Bohm, AUT 

Diophantine Polynomials, D. E. McDougall, Canada 

Graphics World, Currency Change, P. Charland, CAN 

Precise Recurring Decimal Notation, P. Schofield, UK 

Cubics, Quartics - interesting features, T. Koller & J. Bohm 

Logos of Companies as an Inspiration for Math Teaching 

Exciting Surfaces in the FAZ / Pierre Charland 's Graphics Gallery 

BooleanPlots.mth, P. Schofield, UK 

Old traditional examples for a CAS - what's new? J. Bohm, AUT 

What is hiding in Dr. Pest? B. Grabinger, GER 

Truth Tables on the Tl, M. R. Phillips 

Advanced Regression Routines for the TIs, M. R. Phillips 

Where oh Where is IT? (GPS with CAS), C. & P. Leinbach, USA 

Embroidery Patterns, H. Ludwig, GER 

Mandelbrot and Newton with DERIVE, Roman Hasek, CZ 

Snail-shells, Piotr Trebisz, GER 

A Conics-Explorer, J. Bohm, AUT 

and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
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Website for downloading the "Surfer" and other interesting information: 

www . imaginary2 0 0 8.de 

Website with more information about a "Surfer"- competition: 

www . spektrum . de/mathekunst 

Website with results of a former competition (600 graphics): 

www . zeit . de/matheskulptur 

Website for additional information: 

www . algebraic surface . net 



Interruption of calculation for Tl-Nspire CAS software (PC) 

There can be problems cancelling a calculation or interrupting a program with a non-English 
keyboard. You should press the Break (=Pause-) key for some seconds. But it can happen 
that even this does not work. This is a bug and should be resolved with the next release. 


Wolfgang Propper 

... I have tried to generate a family of curves for the G&G-application. It did not work. Can you help 
me? . . . 

I tried (using many tricks including the indirection command) but I failed as Wolfgang. Then I 
wrote to Philippe Fortin and he knew the solution: 


faml\ 

[crx 2 +^,a l { 1,2,3, 4,5,6 }J 
Per tig 

fa ml 3/4 

De fine f a m 1 ( ex, va r, fva ht es) = 

Prgm 
Local i 

For Zjljdimf/vtf/wes) 

expr ( " f 1 1 &str in g ( z )& 11 (x) : = 1 1 &str in g ( ex\ va r= ha In es[ i ] ) ) 
EndFor 
EndPrgm 

i 
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This worked now, thanks to Philippe. And what's about a family of curves given in parameter 
form. Let me try to adapt Philippe's program: 



faml\a m x^+— ,a,{ 1,2, 3, 4, 5, 6} 

terng a 

i 

fampar\ 

[t a-t 2 

Mtt 1 }) 

Fertig 1 

L2 @ 

2/99 

4.08 

0.5 

/ 

fampar 4/5 

De fine f a m p a r( ex, va r, Iva lit es) = 

Prgm 
Local i 

For i, ljdimf/vcj/w^s) 

exprl 11 x " &string(?)& 11 (t) : = 11 &string(ex _ 1 , 1 _ \var -Iva lit es i ) J 
expr( 1 1 y " &str in g ( ? )& 1 1 (t) : = ' ' &strin g ( e x _ 1 , 2 _ \ va r= Iva htesjj) 
c EndFor 

-5.42 

15 5.4 : 

2 EndPrgm 


x3 it 

1 =t 

<t> IE 


y3 {t)=t 2 

0CK6.28 

tstep= 0.13 ^ 

: 


The domain for the parameter must be adjusted in the G&G-application. 


Hubert Langlotz 

. . . Thanks for the rich information in DNL#68. Here is a small program which I used in class. . . 


witrj{l000) 

.02 


Fertig 


wu/y2(iooo) 

.014 

wu?y2(ioooo) 

.0153 

6! 

.015432 

(6) 6 



Jl! 

Simulation of a problem from Probability Theory 

What is the probability to receive numbers from 1 to 6 rollig 6 dies? 

(Bild der Wissenschaft Januar 2008 = Scientific American ) 
wurf2(l0QG)=wurf2(lGG0) 

(Highlight righthand side and evaluate. See examples: 
wurf2(l000)=. 02, wurf2(lQ0Q)=.G2, wurf2(500G)=. 0136, 

In the notes one can only call a function, not a program. wurf2 is the function, wurf is 
the respective program! 


Program and function are among the DNL#69 - files. 
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Ecole de technologie superieure 

Mat 801-01 : Complements de mathematiques (profil Genie Mecanique) Session A-07 

Michel Beaudin 


Devoir #1 

Corrige sur 100, ponderation 15% 

Donne le mercredi 19 septembre 2007 
A remettre le lundi l er octobre au debut du T.P. 

Ce travail peut etre fait en equipe de 2 personnes 

Resolutions symbolique et numerique d’equations/systemes d’equations non lineaires. 
L’importance de pouvoir augmenter la precision. 


1- [25 points] Considerons l’equation algebrique x + yjx + l = 8 . Un logiciel de calcul symbolique a 
donne la solution suivante : 

x = ~— + ^ + 8 ou a = y/: 12V6573+972 et/? = ^12V6573-972 . 

6 6 

Ce qui vaut approximativement 6.07982. 

a) Comment cette reponse a-t-elle pu etre obtenue? 

b) Montrez que la fonction g(x) = 8 — ljx + l satisfait les hypotheses du theoreme du point fixe en 
prenant, par exemple, l’intervalle / = [6, 6.2]. Reconfirmez par la suite la valeur 6.07982. 


Aide : pour a), transformez 1’ equation en une equation polynomiale de degre 3, et utilisez la formule 
de Cardan donnee a la section 0.2 du resume 1 . 


2- [25 points] Considerons l’equationz 26 = 1.01“ z + 3 . 

a) Combien de solutions reelles possede cette equation? Justifiez votre reponse. 

b) Trouvez (par la methode de votre choix) la plus petite solution reelle de cette equation. La 
« methode de votre choix » inclut aussi le solveur numerique de votre systeme de calcul. 

c) Trouvez aussi une solution complexe, en utilisant la methode de Newton a 2 variables, apres avoir 
considere les parties reelle et imaginaire de T equation (une fois la substitution z^x + iy 
effectuee). 
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3- [25 points] L ’ importance de pouvoir augmenter la precision. 


a) 

b) 


20 

Trouvez les 10 solutions reelles de l’equation ]^[(x-£) 
defaut de votre systeme. . . *=i 


— - = 0 . Attention a la precision par 


Soit la fonction f definie par /(x) = x^fx^-Jx + l +yjx — l -2^[x j . Approximez /(lO 6 ) en 
utilisant la precision par defaut de votre systeme. Que trouvez-vous? Augmentez maii^tenant la 
precision et recommencez. Que trouvez-vous? Finalement, montrez que lim/(x) = — — . 

X— >00 4 


2 2 2 2 2 

4- [25 points] Soit la courbe d’ intersection des 2 surfaces x + xyz + y = 3, x +y +2z =4. 
Trouvez les points sur cette courbe qui sont le plus proche et le plus loin de l’origine. 


Aide : notez que « la » courbe d’ intersection des 2 surfaces precedentes est en fait constitute de 2 
courbes fermees (voir figure plus bas). Optimisez le carre de la distance f = x 2 + y 2 +z 2 en 
appliquant la methode des multiplicateurs de Lagrange (voir rappel plus loin). Vous aurez un systeme 
polynomial et, en se restreignant aux solutions reelles, vous devriez trouver 8 solutions (verifiable en 
tenant compte de la base de Grobner associee). 


z 



Quelques rappels en calcul a plusieurs variables 

Rappelons que si w = w(x, y, z) est un champ scalaire a 3 variables, alors son gradient est le vecteur 
dw dw dw 


Vw = 


dx 9 dy 9 8z 


En general, l’equation f(x 9 y,z) = 0 represente une surface dans l’espace et 


le vecteur V/ est perpendiculaire a cette surface. De plus, un vecteur qui est perpendiculaire a la fois 
aux vecteurs v x et v 2 est le vecteur x v 2 , le produit vectoriel de ces 2 vecteurs. 


Lorsqu’on cherche les valeurs extremes de la fonction fix, y, z) sous les contraintes g(x, y, z) = 0 et h(x, 
y, z) = 0, on procede habituellement par multiplicateurs de Lagrange, methode qui dit qu’une condition 
necessaire est donnee par les solutions du systeme de 5 equations a 5 inconnues suivant : 

(*) V/ = AVg + fuVh, g = 0, h = 0 . 
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Solutions du devoir #1, Mat 801-01, session A-07 


Probleme #1 : La reponse a ete obtenue en utilisant la formule de Cardan. En effet: 


1/3 

#1: x + (x + 1) =8 

1/3 

#2: S0LVE(x + (x + 1) = 8, x, Real) 


#3: 


#4: 


1/3 1/3 

(12-76573 - 972) (12-76573 + 972) 

x = + 8 

6 6 

x = 6.079824878 


Voyons d'ou cela est sorti... 

1/3 

#5: (x + (x + 1) = 8) - x 

# 6 : 

1/3 3 

#7: ((x + 1) = 8 - x) 


1/3 

(x + 1) = 8 - x 


#8: 


#9: 


x + 1 = (8 - x) 

3 2 

x = - x +24-X - 192 -x + 511 


3 2 3 2 

#10: (x = - x + 24-x - 192-x + 511) + x - 24-x + 192-x - 511 


# 11 : 


3 2 

x - 24-x + 193 -x - 511 = 0 


Effectuons le changement de variable x=y — 24/3 : 


# 12 : 


24 


y + 


- 24- 


y + 


24 


3 J 


+ 193- 


24 

y + 

3y 


- 511 = 0 


#13: 


y + y + 9 = 0 


#14: S0LVE(y + y + 9 = 0, y, Real) 


1/3 1/3 

(12-76573 - 972) (12-76573 + 972) 


# 15 : 


6 


6 
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Et, puisque x=y + 8, on obtient bien ce que la commande solve du debut a donne. 

1/3 

#16: g := 8 - (x + 1) 

Des graphiques permettent de verifier que g envoie bien fintervalle [6, 6.2] dans lui-meme et que, 
dans ce meme intervalle, la valeur absolue de la derivee de g est majoree par une constante K 
strictement inferieure a 1. Le theoreme du point fixe s'applique done. Prenons fextremite gauche de 
rintervalle: 

#17: FIXED_POINT (g , x, 6) 

#18: [6, 6.087068817, 6.079170206, 6.079884066, 6.079819527, 

6.079825362, 6.079824834, 6.079824882, 6.079824878, 6.079824878, 
6.079824878, 6.079824878] 

Probleme #2: la presence du terme ”+3” fait en sorte que la fonction LambertW de Maple ne peut etre 
utilisee... Une fenetre standard montre 2 intersections (done 2 solutions) mais puisque fexponentielle 
finira par dominer la puissance (pour des valeurs tres tres negatives), une TROISIEME solution reelle 
existe. 


26 -z 

#19: z =1.01 +3 



y ■ 

■10 

■s 

6 

A 


1 1 1 


t 2 i 

X 

1 1 

1 1 1 
A -3 -2 - 

1 

■1 

-2 1 
--4 

1 1 

2 3 4 


La fonction "nsolve" de Derive peut etre guidee avec -oo comme borne inferieure! 

26 -z 

#20: NS0LVE(z = 1.01 + 3, z, -oo, -2) 

4 

#21: z = - 2 . 662524454 • 10 

Une table de valeurs nous aurait permis de constater qu'une troisieme solution existe, situee entre - 
26000 et -27000: 

[ 26 - z 1 

#22: TABLE(Lz , 1.01 + B J , z, -15000, -50000, -1000) 
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4 

- 1.5-10 

108 

3.787675244-10 

64 

6.616171278-10 

4 

- 1.6-10 

109 

2.02824096-10 

69 

1.386693635-10 

4 

- 1.7-10 

109 

9.8100666-10 

73 

2.906392772-10 

4 

- 1.8-10 

110 

4.335958657-10 

77 

6.091553847-10 

4 

- 1.9-10 

111 

1.768453418-10 

82 

1.276738251-10 

4 

- 2-10 

111 

6.7108864-10 

86 

2.675935572-10 

4 

- 2.1-10 

112 

2.386171548-10 

90 

5.608535013-10 

4 

- 2.2-10 

112 

7.998152883-10 

95 

1.175501582-10 

4 

- 2.3-10 

113 

2.540526541-10 

99 

2.463752062-10 

4 

- 2.4-10 

113 

7.682318074-10 

103 

5.163816292-10 

4 

- 2.5-10 

114 

2.220446049-10 

108 

1.082292293-10 

4 

- 2.6-10 

114 

6.15611958-10 

112 

2.268393262-10 

4 

- 2.7-10 

115 

1.642320326-10 

116 

4.754360743-10 

4 

- 2.8-10 

115 

4.227745295-10 

120 

9.964738678-10 

4 

- 2.9-10 

116 

1.052805015-10 

125 

2.088525088-10 


I 4 116 129 

L - 3*10 2.541865828*10 4 . 377372238 *10 J 

26 -z 

#24: f2(z) := z - 1.01 - 3 

#25: RE(f2(x + L- y)) = 0 

f 101 Vx ( C 101 ^ 26 24 2 22 4 

#26: - -COS y-LN +x - 325*x -y + 14950-x *y - 

l 100 J l l 100 )) 

20 6 18 8 16 10 
230230 • x *y + 1562275-x *y - 5311735-x *y + 


14 12 12 14 10 16 

9657700 • x *y - 9657700-x *y + 5311735-x *y 


8 18 6 20 4 22 2 24 26 

1562275-x -y + 230230-x -y - 14950-x -y + 325-x -y - y -3=0 


#27: IM(f2(x + L- y)) = 0 
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#28: 


101 


100 


Vx ( ( 101 ^ 

•SIN y • LN 

J l l 100 )) 


24 22 2 

+ 2-x-y-(13-x - 1300 • x -y + 


20 4 18 6 16 8 14 10 

32890 • x .y - 328900-x -y + 1562275-x -y - 3863080-x -y 


12 12 10 14 8 16 

+ 5200300 • x -y - 3863080-x -y + 1562275-x -y 


6 18 4 20 2 22 24 

328900 • x -y + 32890-x -y - 1300-x -y + 13-y ) = 0 


En tragant les 2 courbes (#26 et #28) dans une meme fenetre, une infinite ^intersections (solutions 
complexes de Fequation fl(z) = 0) se presentent. Nous allons choisir le couple (0.12, 1) 



#29: NEWT0NS([RE(f2(x + L- y)), IM(f2(x + L-y))], [x, y] , [0.12, 1], 10) 


0.12 

1 

0.1356956448 

1.088548039 

0.1311422562 

1.061873117 

0.1281462587 

1.049333933 

0.1272824115 

1.04710648 

0.1272412487 

1.047049788 

0.127241196 

1.047049764 

0.127241196 

1.047049764 

0.127241196 

1.047049764 

0.127241196 

1.047049764 

. 0.127241196 

1.047049764 . 


On aurait pu utiliser Newton a une seule variable complexe: 

#31: NEWT0N(f2(z) , z, 0.12 + L, 10) 

#32: [0.12 + L , 0.1356956449 + 1. 08854804- i , 0.1311422563 + 1 . 061873118 • L , 

0.1281462588 + 1. 049333934- L , 0.1272824116 + 1 . 04710648 • L , 
0.1272412488 + 1 . 047049788 • i , 0.127241196 + 1. 047049764- i , 
0.127241196 + 1. 047049764- i , 0.127241196 + 1 . 047049764- L , 
0.127241196 + 1 . 047049764- i , 0.127241196 + 1. 047049764- l] 
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Probleme #3: pour a), "nsolve" ou "nsolutions" est rapide mais des solutions entieres sont 
evidemment impossible... 


#33: expl 


(20 

n (x - k) 

U=i 


19 

x 


7 


10 


#34: NSOLUTIONS (expl , x, Real) 

#35: [1, 2, 3, 4, 5, 20.78880561, 6.999745978, 6.000005824, 

8.006075438, 8.928803402] 

#36: APPROX(NSOLUTIONS(expl, x, Real), 20) 

#37: [1, 2, 3, 20.78880561, 8.928803402, 8.006075438, 6.999745978, 

4.000000000, 6.000005824, 4.999999939] 

#38: APPROX(NSOLUTIONS(expl, x, Real), 30) 

#39: [1, 2, 20.78880561, 8.928803402, 4.999999939, 8.006075438, 

6.000005824, 2.999999999, 6.999745978, 4.000000000] 

#40: APPROX(NSOLUTIONS(expl, x, Real), 40) 

#41: [1, 8.928803402, 20.78880561, 2.000000000, 2.999999999, 

8.006075438, 4.000000000, 6.000005824, 6.999745978, 4.999999939] 
#42: APPROX(NSOLUTIONS(expl, x, Real), 50) 

#43: [6.000005824, 8.006075438, 6.999745978, 4.999999939, 20.78880561, 

2.999999999, 0.9999999999, 8.928803402, 2.000000000, 4.000000000] 
#44: APPROX(NSOLUTIONS(expl, x, Real), 60) 

#45: [4.000000000, 6.999745978, 2.999999999, 0.9999999999, 8.006075438, 

6.000005824, 20.78880561, 8.928803402, 2.000000000, 4.999999939] 
#46: exp2(x) := x-Jx'(J(x + 1) + 7(x - 1) - 2-,/x) 


Pour calculer la limite de exp2(x) lorsque x tend vers l'infini, on peut, par exemple, remplacer x par 
l/h, simplifier et laisser h tendre vers 0 par la droite. La regie de l'Hospital s'applique alors et, en 2 
coups, on trouve la valeur -1/4 (voir #54 plus bas). Derive trouve cette valeur sans nous indiquer 
comment... 


#47: exp2(x) 

lim exp2(x) 
#48 : x-co 


1 


#49: 


4 


12 
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#50: exp2 


h J 


7 


1 - h 


7 


h + 1 


■ SIGN(h) 


#51: 


3/2 


#52: h_ :e Real (0, ~) 


#53: 


7 

' 1 - h_ y 

+ 7 

' h_ + 1 V 



V 

< h_ > 


s. h_ ) j 


#54: 


■ SIGN(h_) 


3/2 2 

h_ h_ 

V (1 — b_) + 7(h_ + 1) - 2 


#55: 

#56: 


exp2(10 ) 


#57: 

#58: 


APPROX(exp2(10 ), 20) 


-0.2500000003 


#59: 

#60: 


APPROX(exp2(10 ), 50) 


-0.2500000000 


Probleme #4: notez que le graphique a ete produit en utilisant RK pour le systeme d'E.D. provenant 
du produit vectoriel des 2 gradients... 


.: Ldi: 


2 2 2 2 2 2 
#61: Ldist2 := x + y + z , contl :=x + x-y-z + y - 3, cont2 := x + 


2 2 1 
y + 2-z - 4 J 

Les 5 equations provenant de la methode des multiplicateurs de Lagrange sont rapidement obtenues 
par Derive: 

#62: APPEND_C0LUMNS(GRAD(dist2) - A-GRAD(contl) - p •GRAD(cont2) , 

[contl], [cont2]) 

[ - x-(2-A + 2-p - 2) - A-y-z, - A-x-z - y-(2-A + 2-p - 2), z-(2 - 

2 2 2 2 2 1 
4*p) - A*x*y, x + x-y-z + y -3, x +y + 2*z - 4j 


# 63 : 
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#69: [1.732050807, 1.732050807, 1.732050807, 1.732050807, 1.963186685, 

1.963186685, 1.963186685, 1.963186685] 


Done, les points (1, 1, 1), (1, -1, -1), (-1, 1, -1) et (-1, -1, 1) sont tous a une distance de 73 = 1.732 
unite de l'origine (et sont done les plus proches) tandis que les 4 autres points sont tous a une distance 
de 1.96 unite de l'origine et sont les plus eloignes. 



Two examples of pedalcurves 
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An Optimization Problem: From a Spreadsheet to a CAS 

Josef Bohm, Wiirmla, Austria 


In the MNU Journal 60/6 from 2007 (MNU = Mathematisch-Naturwissenschaftlicher Unterricht = 
Education in Maths and Science ) Dr. Helmut Brunner presented an optimization problem and its solu- 
tion using the Solver from MS Excel recommending the spreadsheet as the appropriate tool for solving 
this and similar problems. 

I read his contribution with interest - one reason is, that I have known Dr. Brunner since many years 
very well from my times as teacher. He was headmaster of a College of Business Administration in 
Modling, Lower Austria, for many years, speaker in numerous seminars, textbook author and well 
known in the math community. I appreciate him as an excellent mathematician. He started very early 
using the programmable pocket computers for math teaching and he was very successful as head of the 
teachers working group in our federal state propagating the importance of probability theory and sta- 
tistics in math education. 

But I cannot agree with Dr. Brunner (Helmut) praising a spreadsheet as “the“ appropriate tool for tack- 
ling this problem. 

The problem: Find the pairs of points P(xi,yi) and Q(x 2 ,y 2 ) having minimal or maximal 
distance from each other with P being a point on the ellipse x 2 + 4y 2 = 4 
and Q being a point on the circle (x - 2) 2 + (y - 2) 2 = 1 . 


This extremal value problem with four variables and two side conditions will be solved by using 
Lagrange multipliers: the extended main condition must be derived partially wrt to the four variables 
xi, y u x 2 und y 2 . These derivatives will be equated 0 and form together with the side conditions a 
nonlinear system of equations with 6 unknowns. 


Main condition: 



Side condition 1 x 2 + 4 y 2 -4 = 0 


(eql) 


Side condition 2 


(x 2 -2) 2 +(y 2 -2) 2 -\ = 0 


(eq2) 


Helmut takes the square of the distance which gives the following main condition extended by the 
multipliers: 

mcl = (x 2 - x 1 ) 2 + (y 2 - y 1 ) 2 + A • (^x 2 + 4y 2 -4j + ju- |(x 2 - 2) 2 + (y 2 - 2) 2 - 1 j 

The partial derivatives give equations (3) to (6): 

u mcl c mcl 

= -2 (x 2 - jfj ) + 2Zx l = 0 (eq 3) = -2 ( y 2 - y l ) + 8 Zy l = 0 (eq 4) 

c.v, dy x ' ' ' 

o mcl (j mcl 

— — = 2 ( x 2 -x l ) + 2/i{x 2 -2) = 0 (eq 5) — = 2(y 2 - y x ) + 2/^{y 2 -2) = 0 (eq6) 

GX 2 2 
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Helmut eliminates (by hand) the multipliers X and //. (Could be a nice task for students to repeat basic 
skills in manipulating expressions!) The remaining system of equations is complicated enough to resist 
all efforts to solve it by hand. 


Xj 2 + 4y^ = 4 

(eq 1) 

(x 2 - 2) 2 + (v 2 - 2) 2 = 1 

(eq 2) 

4 v, (x 2 - x, ) = x, ( v 2 - v, ) 

(eq 3) 

( y 2 - 2 ) 0*2 - *0= 0*2 - 2 ){yi - Vi ) 

(eq 4) 


Helmut does not even try to solve this system but he uses Excel and its Solver and finally obtains after 
setting appropriate initial guesses (which ones??) for the - hidden in the dark - iteration algorithm in a 
straight forward way two solutions. I show the mimimum solution: 



One can read off the coordinates of the points 
P and Q and the minimum distance. For finding 
the maximum one has to change from Min to 
Max and change also the initial values. 


A4 

£ =(l 

(C2-A2) A 2+(D2-B2) A 2) A (1 /2) 

A 

B 

C 

D 

L 

1 

XI 

Y1 

X2 

Y2 


2 

1,33562112 

0,72111969 

1 ,56693693 

1,09361245 


3 






4 

0,41880109 1 

1 




5 

4,00000035 ' 

4 




6 

0,99999979 

1 





As you can see, he does not solve the system, but tries to find the minimum/maximum under the given 
side conditions. You can also try to solve the complete - reduced - system using the solver: 


A 

B 

C 

D 

XI 

Y1 

X2 

Y2 

1,33564120926219 

0,721110075022646 

1 ,56693736779241 

1 ,09361204343036 





=((C2-A2) A 2+(D2-B2) A 2) A (1 /2) 




=A2 A 2+4*B2 A 2 

4 



= (C2-2) A 2 +(D2-2) A 2 

1 



=4*B2*(C2-A2)-A2*(D2-B2) 

0 



= (D2-2)*(C2-A2)- (C2-2)*(D2- B2) 

0 




(The Excel files are among the downloadable files.) 
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Doubtless the Solver is an excellent numerical tool and I have used it very often - when working with 
Excel. But nowadays in many schools other - mathematical - tools are at disposal. First of all I have 
CAS in mind - what else? Most of the Computeralgebra Systems have built in the Groebner Bases 
(Learn more about Groebner Bases in DNL#68!) and are able to solve polynomial equation systems. 
Application of these systems set the main activity of the students into the mathematical definition and 
description of the problem and then into the interpretation of the solution (if there is one). Of course, 
the students need a certain familiarity with ’’their” CAS. But when students have to solve problems 
like this they will have acquired some experience. 

In the following Ell show the complete procedure using DERIVE (followed by the solution process 
performed on the handheld and finally with TI -Nspire). 


DERIVE allows parallel presentation in form of „abstract formulae 46 and „visible geometric objects 44 . 
This is an enormous advantage for teacher and student. 


#1: InputMode := Word 

2 2 

#2: x + 4-y = 4 

2 2 

#3: (x - 2) + (y - 2) = 1 



We set up the main condition (taking the square of the distance) and the two side conditions (the 
points are lying on the two respective curves) and introduce the Lagrange multipliers. Having done 
this, the mathematical problem is solved. We can leave the solution of the system of equations to our 
tool, a CAS. I try to keep the problem as general as possible and introduce parameters m, n and r for 
the coordinates of the center of the circle and its radius. (This proves to be helpful for further investi- 
gations.) 

2 2 


#4: 

me := (x2 - 

- xl) + (y2 - yl) 

#S : 

2 

scl := xl 

2 

+ 4-yl -4=0 

#6: 

sc2(m, n, 

2 2 2 
r) := (x2 - m) + (y2 - n) - r = 0 

#7: 

mclfm, n, 

r) := me + X-LHS(sel) + p ■ LH5(sc2(m ! 


The MNU-Problem: m = n = 2, r = 1 . 

#8 is for the plot and #9 gives the solutions. 

#8: [scl, sc2(2, 2, 1)] 


#9: 


f d d d d 

SOLUTIONS mcl (2 , 2 , 1) a mcl(2, 2, 1) a mcl(2, 2 , 1) a mcl (2 , 

Id xl d x2 d yl d y2 


COL [1, 2, 3, 4] 


2, 1) a scl a sc2(2 , 2, 1) 


[xl, yl, x2 , y2 , X, p] 
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Expression #9 contains the complete system. As we are not interested in the values for X and // output 
of the first four columns will be sufficient: 


1. 

3 SB 64 

0. 

72111 

2.43301 

2 

.90136 

-1. 

92051 

-0. 

279113 

2 .36453 

2 

.50257 

1. 

3 SB 64 

0. 

72111 

1.5669E 

1 

.09661 

-1. 

920B1 

-0. 

279113 

1.13546 

1 

.49742 

3.65356 

# 4.06302 ■ i 

-2.16447 

+ 1.80843 -L 

3.65358 # 4.06302 ■ l 

-2.16447 

+ 1.60643 -L 

3.65356 

- 4.06302 ■ L 

-2.16447 

- 1.S0S43 ■ L 

3.85356 - 4.06302 ■ 1 

-2.16447 

- 1.80843 -L 

2.9341 + 1.44216 - L 

-0.SS7665 

+ 1.19173 -L 

2.16497 + 0. 62 4219 -l 

0.629455 ■ 

#■ 0. 0679741 -i 

2.9341 - 

■ 1.44216 - L 

-0.SS7665 

- 1.19173 -l 

2.16497 - 0. 62 4219 -L 

0.629455 ■ 

- 0 .0879741 - L 

2.9341 hj 

■ 1.44216 - t 

-0.SS7665 

+ 1.19173 -L 

1.33502 - 0 .624219 - L 

3.17054 ■ 

- 0. 0879741 - L 

2.9341 - 

■ 1.44216- l 

-0.SS7665 

- 1.19173 -L 

1.33502 + 0 .624219- l 

3.17054 ■ 

i- 0. 0879741 - L 

1.47974 i 

■ 0. 729097 -l 

0 .S3 114 - 

0 .324517 ■ L 

1.47974 + 0. 729097 -l 

0.63114 ■ 

- 0 .3245 17 ■ t 

1.47974 - 

■ 0. 729097 -L 

0 .33114 + 

0 .324317 ■ L 

1.47974 - 0.729097- l 

0.63114 ■ 

i- 0 .3245 17 ■ l 


We don't need the complex solution, so we include the restriction on real values only. As I am plan- 
ning to perform this procedure not only once, I define a function to deliver the real solutions for vari- 
ous parameters m, n and r. Then I apply this function for the MNU-paper parameters: 


#11: sols(m, n, r) := SOLUTIONS 


4 


i 


d d d 

mcl(m, n, r) a mcl(m, n, r) a me 1 ( m , n, 


fd xl 


d x2 


d yl 


r) a 


d y2 


mcl(m, n, r) a scl a sc2(m, n, r), [xl, yl, x2, y2, X, p], Real COL 


)) 


[1, 2, 3, 4] 


sols(2 P 2. 

1) 



' 1.36564 

0.72111 

2.43301 

2.90138 

-1.92051 

-0.279113 

2.86453 

2.50257 

1.38564 

0.72111 

1.56698 

1.09861 

. -1.92051 

-0.279113 

1.13546 

1.49742 


Here we find the second important difference to the Excel-Solver-solution. (The first one is - in my 
opinion - the representation of the equations in DERIVE which is very close to the hand written pro- 
cedure. I really appreciate the concept of the relations between cells in a spreadsheet and all the other 
valuable benefits of a spreadsheet, but this is a better job for a CAS.) Four solutions are presented in 
one single step and we need not enter any guesses. The criteria for determining the type of the extre- 
mal values are so difficult, that they are even rarely given explicit in textbooks. It is very easy to check 
the validity of the solutions by visualisation and inspection. We plot the solution points together with 
their connecting segments and we (let) calculate the respective distances. This can be done without 
any problems. 

The third - important difference - is the easiness of visualisation of the solutions which can inspire for 
further investigations. 
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1- Solution 

1.38564 0.72111 

2.433 2.90138 


#14: 


-/CSUBSTfmCp [xl P yl p x2 ? y2] r (sol s(2 r 2, 1)) )) 
#15 : 1 


#16: 2.418801041 

2. Solution 

-1.92227 -0.27915 

2.86491 2.50258 J 


#17: 


# 20 : 


-/(5UB5T(mc . [xl ? yl, x2 p y2] p (sols(2 f 2 f 1)) )) 
#21: 3 


#22: 0.4188010415 


The Minimum 


4. Solution 



#23: 

" -1.92227 
_ 1.13547 

-0.27915 ' 
1.49742 _ 


./(5UBST(mc , 

#24: 

[xl ( yl, x2, y2 ] , Csols(2 r 2, 1)) )) 

4 

#25: : 

3.534843535 




-/(5UB5T(mc ? [xl P yl P x2 f y2] P Csolst2 r 2 r 1)) )) 
#18: 2 

#19: 5.53484358 

The Maximum 

3. Solution 

1.38564 0.721109 

1.56698 1.09861 





An experienced DERIVE user will find the values for all distances belonging to the pairs of solution 
points in one single step: 


#26: 


VECTOR 


( ( 

Cv 
l l 3 






sols(2 p 



#27: [2. 418801041 P 5. 5 348435 8 , 0 .4188010415 p 3.534843535] 

The 2 nd pair leads to the maximum and the 3 rd pair to the minimum. 
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The CAS-solution offers an interesting question for the students: Where do the other solutions come 
from? Do they have any interpretation? (I will come back to this question later in this paper.) 

As you might know the algorithm for solving systems of polynomial equations (Groebner bases) is 
also implemented in symbolic pocket calculators as TI-92+, Voyage 200 and Tl-NspireC AS. The next 
figures show the procedure performed using Tl-NspireCAS. Capacity of the devices is restricted by 
their smaller memory, but often one can help himself. 




Even the reduced system from MNU6/60 proves to be too complex for the TI-92+/Voyage 200 and 
you will receive the message about too little memory. With little support, which addresses the stu- 
dents' competence in handling the device one can ’’persuade” even the TI-92+ to release the solution 
what shall be demonstrated in the following: 
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I have to use variables a, t , u and v, because yl andy2 are system variables on the TI-92 and Voyage 
200. The system consisting of the four equations from above can be reduced once more leading to a - 
nonlinear - system with three unknowns. I substitute for u and this system can now be solved com- 
pletely. The device gives 6 triples of solutions: 

s = 2,t = v = Q 
s = -2, t = v = 0 

s — 1,385641, t = 0,721 119, v = 1,098612 
s = 1,385641, t = 0,721 1 19, v = 2,901388 
s = -1,920516, t = -0,2791 13, v = 2,502578 
s = -1,920516, t = -0,2791 13, v = 1,497422 

The first two solutions don't make sense because of t = 0, the other four solutions are the same as in 
the DERIVE calculation. The missing value for u can be found by substitution. Then it is easy to cal- 
culate the distances. Maximum and minimum are found by inspection. Plotting is not so easy as it is 
using DERIVE. 

In addition to the discussion about the ’’other” solutions hopefully the interest to inform about Grobner 
Bases and the Buchberger algorithm will be roused. So a ’’Black Box” can be coloured - not bright 
white - but at least grey. 

Dr. Brunner's method is justified when students only have access to and experience with Excel. But it 
is my strong opinion that if a CAS can be used this should be preferred. 


This is what - except the remarks on further investigations and some smaller changes - 1 wrote in my 
comments to Dr. Brunner's article which appeared - without showing the TI-92 & TI -Nspire treat- 
ments - in MNU 61/1. As I am here not restricted in space in the DNL I can proceed . . . 
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My first intention was to solve the problem without using the multipliers because they are not part of 
the curriculum in secondary schools. Students might be inspired by the plots to try the following ap- 
proach: 

Find pairs of points on the two curves with parallel tangents and with the connecting segment perpen- 
dicular to the segments. In this case we need implicit differentiation of the curves. 


x 2 + Ay 2 = 4 


(x-2) 2 +(j-2) 2 =1 


2x + 8 yy' = 0 — » y' = - 

2(x-2) + 2(y-2)y' = 0^>y' = - 



2 

~2 


Points P(xi,yi) and Q(x 2 ,y 2 ) are lying on the ellipse and circle, the slopes in P and Q must be equal and 
the slope of segment PQ is the negative reciprocal value of the slope of the tangents. This leads to four 
equations (without needing the multipliers because it is no extremal value problem any longer!): 


x 2 + 4y 2 -4 = 0 (eql) 

(x 2 -2) 2 +(y 2 -2) 2 -1 = 0 (eq 2) 


x 2 -2 

y 2 ~ 2 


4 Ji 


4 y, 0 


2 ) = x i(y 2 - 2 ) 


(eq 3) 


y 2 -y\ 

X 2 -Xj 


y 2 zl 

x 2 -2 



2 )(^2 - J ^ l ) = (>"2 - 2 )(^2 -^ l ) 


(eq 4) 


Without looking back in the article we leave the solution again to the CAS: 


#26: [eql := scl, eq2 := sc2(2 r 2 r 1)] 

#29: eq3 4-yl-(x2 - 2) = xl-(y2 - 2) 

#30: eq4 := (x2 - 2)-(y2 - yl) = (y2 - 2)-fx2 - xl) 

#31: S0LUTI0NS(eql a eq2 a eq3 a eq4, [xl, yl, x2 , y2] , Real) 


r -1,92051 

-0.279113 

1.13546 

1.49742 

1.38564 

0.721110 

1.56698 

1.09861 

-1.92051 

-0.279113 

2.86453 

2.50257 

_ 1.38564 

0.721110 

2.43301 

2.90138 


Compared with the solution from above we don't find any difference - 
- except the order of the solutions: 


' 1.38564 

0.72111 

2.43301 

2.90138 ‘ 

-1.92051 

-0.279113 

2.86453 

2.50257 

1.38564 

0.72111 

1.56698 

1.09861 

_ -1.92051 

-0.279113 

1.13546 

1.49742 . 


Comparing the system with the reduced system from page 1 6 we find that both systems are nearly the 
same. Only in eq 3 the variables x\ and yi are both substituted by 2. Any idea why? 
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Before talking about this we will use our generalized functions to investigate other cases. What will 
happen if the circle lies completely within the ellipse, what if both curves will intersect, . . . 

We choose m = -3, n = 2 and r = 7: 


#39: 

#40: 

#41: 

#42: 


Notation Dig its := 5 
[scl p sc2(-3 P 2 r 7}] 
sols (-3 r 2, 7) 


' 1.9516 

-0.2186 

3.388 

-0.86222 

-1.7254 

0.5057 

-7.5 427 

7.3257 

-1.7254 

0.5057 

1.5427 

-3.3257 

. 1.9516 

-0.2186 

-9.388 

4.8622 


#43: 


1.9516 -0.2186 

3.388 -0.86222 


-1.7254 0.5057 

-7.5427 7.3257 


-1.7254 0 .5057 

1.5427 -3.3257 


[vf(y - V ) t- (v - v ) ] F v ? sols(-3 f 2 r 7)1 
L L 3 1 4 2 J j 


#44: VECTOR 

#45: [1.574, 8.964, 5.0359, 12.425] 

#46: 


T 1.9516 -0.2186 ' 


" 1.9516 -0.2186 T 

_L 3.388 -0.86222 . 

? 

. -9.388 4.8622 J. 


1.9516 -0.2186 11 

-9.388 4.8622 J. 


There are again 4 solutions. 

Plot of #43 with medium sized points and dotted 
connection lines show all solutions, #46 with 
large sized points and solid connections give the 
presentation of maximum and minimum. 

Any idea to reduce the system of equations once 
more? 

In the next experiment we have intersecting 
curves: 



m = -3, n = -2. r = 4 

#47: [sclp sc2(-3, -2, 4)] 

#48: sols (-3 f -2, 4) 


-1.7254 

-0.5057 

-5 .5958 

-5.0432 ' 

-0.33802 

0.98561 

-0.33802 

0.98561 

1.9516 

0.2186 

-6.6503 

-3.6355 

-1.7254 

-0.5057 

-0.40416 

1.0432 

1.9516 

0.2186 

0.65033 

-0.36444 


L 0.84742 


■0.90579 


0.84742 -0.90579 
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#51: VECTOR 


f f 

(v 
l l 3 


V ) + (v - V } l r Vp sol!s(-3 , -2, 4)1 

1 4 2 J J 


T 1.9516 0.2186 ' 


" 1.9516 

0.2186 1" 

_L -6.6503 -3.6355 . 

r 

. 0.65033 

-0.36444 J. 


#52: [5.964, 0 P 9.4259, 2.0359, 1.4259, 0] 


#53: 


This problem leads to six solutions of the system. 

Two of them show the final function value 0. 
What is the interpretation of them? 


Any idea to reduce the system of equations once 
more? 


What if the circle lies completely inside of the 
ellipse? 


m = n = 1/4 P r = 1/2 

r i 



#54: 


scl, sc2 


#55: sols 


#56: 


1.9899 


-■ -■ -)1 
4 4 2 JJ 

-■ -1 
.4 4 2 J 


#58: VECTOR 


1.9899 

0.3074 

0.36419 

0.3074 

( ( 

I c 


-0.1001 

0.74017 

0.15137 ' 

-0.98328 

0.20389 

0.74787 

-0.071402 

0.74495 

0.32086 

-0.071402 

-0.24495 

0.17913 

-0.1001 

-0.24017 

0.34862 

0.98811 

0.21122 

-0.24849 

-0.98328 

0.2961 

-0.24787 

0.98811 

0.28877 

0.74849 . 

2 

- v ) + 

3 1 

(v - v ) 

4 2 . 

j. v, sols 


#59: [1.2748, 1.7385, 2.7677, 1.7677, 2.2748, 1.2403, 0.73855, 0.24034] 


# 60 : 


T 0.3074 0.98811 " 


[ -1.9948 -0.071402 11 

_L 0.28877 0.74849 . 

T 

. 0.74495 0.32086 J. 



This gives us eight solutions. Is it possible to have more than eight solutions? The plot demonstrates 
very clear what we sould have known for a while: 


The connecting lines of all pairs P and Q are passing the center of the circle (which is a consequence 
of the fact that the segment is perpendicular to the parallel tangents). This is presented algebraically by 
substitution of x 2 andy 2 by 2 in the above equations on pages 16 and 22. 


We can again redesign the solution process (which leads to a reduction the system of the equations): 

Find the point(s) P on the ellipse where the tangent is perpendicular to the line PM (M = center of the 
circle). 

x 2 +4y 2 -4=0 (eql) 

^- = — -> x l (y l -n) = 4y l (x l -m) (eq 5) 

x, — m x, 


#61: eqS(m P n) xl-(yl - n) = 4-yl-(xl - m) 

#62: sols2(in P n) := SOLUTIONS (eql a eqS(m r n) r [xl r yl] , Real) 


#64: 


# 66 : 


sols2(2 . 2) 


" -1.92051 

-0.279113 

l 1.35564 

0.721110 

f 1 

sols2 P 

L 4 

t) 

" 1.95995 

-0.100101 

-1.99469 

-0.0714024 

0.307407 

0.965116 

. 0.364199 

-0.963250 


We have to find the intersection points of all lines PM with the ellipse in order to find the complete 
solutions and then select the points giving maximum and minimum distance. 

I hope that you enjoyed my CAS-treatment and the extensions of the given Excel-Solver-application 
and that you share my preference of using CAS - symbolic calculation and plotting simultaneously. 
Using only the spreadsheet we miss a lot of opportunities to do mathematics. 
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Analysis of Variance - ANOVA 


MacDonald R. Phillips 
don.phillips@gmail.com 


#1: Notation := Decimal 

#2: Precision := Approximate 

In general, the statistical technique analysis of variance, ANOVA, is used to test the hy- 
pothesis that the means of three or more populations are equal. The routines below can be 
used for one-way analysis of variance, ANOVA, two-way analysis of variance, 
TwoWayANOVA, and two-factor analysis of variance with interaction, TwoFactorANOVA. 

1 . One-Way Analysis of Variance (ANOVA) 

ANOVA takes a matrix of three or more rows, each row representing a sample from a differ- 
ent population. The samples do not have to be of equal size. ANOVA tests whether the 
means of the samples are equal. 


Example 1: Quality-awareness examinations were given to 6 employees at three different 
manufacturing plants. Based on the scores below, are the average scores the same in each 
plant? 




' 85 

75 

82 

76 

71 

85 ' 

#3: 

data := 

71 

75 

73 

74 

69 

82 



. 59 

64 

62 

69 

75 

67 . 


#4: ANOVA(data) 

One-Way Analysis of Variance 
f Ho: Means Equal 

Ha: Means Not Equal 

<x: 0.05 

L F(<x): 3.682320331 


Source 

DF 

ss 

MS 

F 

P(F) 

Between 

2 

516 

258 

9 

0.002702899474 

Wi thi n 

15 

430 

28.66666666 



. Total 

17 

946 





Conclusion: Reject Ho 
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The ANOVA output shows Ho, Ha, a (the level of confidence) and the F value of a. Next it 
shows the ANOVA table with DF (the degrees for freedom), SS (sum-of-squares), MS (mean 
sum of squares), the F value for the test statistic, and P(F) (the probability of F. Since P(F) is 
less than a, the conclusion is to reject Ho and accept the Ha that not all the means are equal. 
This means that some of the means may be equal, but not all of them are. Other tests are 
needed to determine if all the means are not equal, or if some are, which ones they are. 

The test may also be plotted showing the curve of the F-distribution, the rejection region, and 
the value of the test statistic. 

#6: ANOVAPlot(data) 



How to do it without Derive? 

The Tl-Voyage and Nspire offer a bundle of statistics tools with ANOVA among them. You 
can follow how to treat the first example using the Stats/List Editor-Application of the Voyage 
200/ TI-89/92. 

Please compare the results. Plotting the F-distribution and the rejection area is not so easy, 
but it can be done using the built-in fpdf-function together with the shade-option in the Graph- 
Window: 


r FIT Y Y F]:t Y FHt 

Tools Plots List Calc 

f F ST wxa/ F FT Y 

□ i str IC5H I nts 1 

ex 1_1 

ex 1_2 

ex 1_3 

ex 

1 : Z-Test... 

2: T-Test... 

3: 2-SanpZTest... 

4: 2-SampTTest... 

5: 1-PropZTest... 
6:2-PropZTest... 

7: Chi 2 GOF... 
8:Chi2 2-way... 

9: 2-5ampFTest... 
fl: LinRegTTest 
B: MultRegTests... 1 

35 

71 

59 

75 

75 

64 

32 

75 

62 

76 

74 

69 





exl_3[3]=62l 


MAIN RAD EXACT FUNC 3/ 9 


r F1 t Y FEt Y F3t Y FHt - 

|Tools|Plots|List|Calc 

f FSt Y FbT y F ?tY 1 

|D i str |Tests 1 1 nts | 

ex 1_1 

ex 1_2 

ex 1_3 1 







35 , 

75 


Choose Input Method 



32 

76 

71 

35 \ 

Data Input Method: Data-* 
Number of Groups: y* 


CEnter=0K h <TESC=CANCEO 





r 






exl_3[7]= 

USE * AND + TD OPEN CHOICES 


\ FIT Y Y f:=t T FHt T Y FbT Y Y 'T 

|Too 1 s |P 1 ois |l i si |c a 1 o |P i sir |Tests 1 1 nts | | 

ex 1-1 I ex 1_2 lexl_3 I -| 


35 

75 
32 

76 
71 
35 


exl_3E7] = 


f Analysis of Variance 


List 1: 

ex 1_1 


List 2: 

exl_2 


List 3: 

exl_3l 


, CEnter=0K } CESC=CRNCEL^ ; 


f fit 

|Tool 

f Analysis of Variance "> 

pn 

F -9.GGGGGG 

P Ualue =.0G27G3 

FACTOR: 

df =2.000000 

S3 =516.000000 

MS =25S.GGGGGG 

ERROR: 

4-df =15.GGGGGG 

^£Enter=0K 5 

exl_ 

pi ist 

35 

75 
32 

76 
71 
35 

3.659 

3.659 

upl 


MAIN RAD EXACT FUNC fi/ b 
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iToo ll 

exl. 

35 

75 
32 

76 
71 
35 


upl 


f 

Analysis of Variance 


■rdf 

=2.000000 


ss 

=516.000000 


MS 

=253.000000 


ERROR: 

df 

=15.000000 


SS 

=430.000000 


MS 

=23.666667 


Sxp 

^£Enter= 

=5.354126 
=0K j 

J 


D 

pi ist 
3.653 
3.653 


r FIT Y F£t Y FjT Y FHt " 

|Tools|Plots|l_ist|Calc 

f F Et Y FfiT Yf?tY } 

|0isir|Tesis|lnts| 

exl_2 

exl_3 

xbarl... 

lowl i... 

upl ist 


71 

75 

73 

74 
63 
32 

53 

64 

62 

63 

75 

67 

73.000 

74.000 

66.000 

74.341 

63.341 

61.341 

33.653 

73.653 

70.653 





xbarlist=<!79 . , 74 . ,66 _> 

MAIN RAD EXACT FUNC H/ b 


The Numeric Solver provides the bound of the rejection area. Then we can plot the distribu- 
tion function and shade the rejection area. 



fTi-^mnY Y fs yfh Yf^ty Ffi-r y *:* y ) 

| t |Zoon |Ed i 1 1 |fl 1 1 |sty 1 e |r'i .. | 

J-LOTl' 

Plot 2:LH 
Plot 1 
■^yl=tistat 
y2=l 

y 5 = 

y4= 

y5= 

y6= 

y£= 

y§= 

□ >: : ci :-':cb 

■ X-Cl v: CE 

■ fpdf(x,2. 

15) 

w2<x>= 

STAT 

RAD AUTD 

FUNC 








STAT RAD AUTD FUNC 


I adapt the Window-settings to obtain a better visualisation of the rejection area. The shading 
needs the lower bound 3.68. 


fr f^Izoonjll racelRegraphl 




1 : Ualue 
2: Zero 
3: Mininun 
4: Maxinun 
5: Intersection 
6: Derivatives ► 
7: Yf <x)dx 
3: Inflection 
3: Distance 
fl: Tangent 
E: Arc 











STAT 

RAD AUTD 


FUNC 



There is a second example illustrating the One-Way-Analysis followed by the Tl-Nspire solu- 
tion. 

Example 2: 


#17: example 


9 

ii 

7 

8 

9 

9 

10 ■ 

11 

12 

8 

9 

10 

11 

10 

12 

IB 

11 

9 

10 

12 

12 

15 

17 

15 

10 

16 

14 

12 . 
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#18: ANOVA(exampl e) 


#19: 


One-Way Analysis of Variance 


Ho: 

Means Equal 

Ha: 

Means Not Equal 

or. 

0.05 

F(<x) : 

3.008786556 


Source 

DF 

SS 

MS 

F 

P(F) 

Between 

3 

102.2857142 

34.09523809 

12.18723404 

-5 

4.813001956*10 

Wi thi n 

24 

67.14285714 

2.797619047 



Total 

27 

169.4285714 





Conclusion: Reject Ho 


#20 : AN0VAP1 ot (exampl e) 



/7: = {9,ll,7,8,9,9,lo} 

{9,11,7,8,9,9,10} 

A 

/2: = { 11,12,8,9, 10,11,10} 

{ll,12,8,9,10,ll,10} 

/3:={ 12,13,11,9,10,12,12} 

{12,13,11,9,10,12,12} 

!4:={ 15,17,15,10,16,14,12} 

{15,17,15,10,16,14,12} 

ANOVA 11,12,13,14 , 0 


Per tig 

j slat. results 

"Titer 

"ANOVA" 



11 F" 

12.1872 



"PVal" 

.000048 



"df" 

3. 



11 SS" 

102.286 



"MS" 

34.0952 



"dfError" 

24. 



"SSEiror" 

67.1429 



"MSError" 

2.79762 



"sp" 

1.67261 



"CLowerList" 

"{-r 



"CUpperList" 




"xList" 



invF(.95,3,24) 

| 


3.00879 " 
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Compare again the results. Working with Tl-Nspire is very similar with working on the TI-92 
or Tl-Voyage 200. The CAS-machine which does the work is the same one. 

The additional advantage of Don's tool is that the answer if the hypothesis should be re- 
jected or not is given in the result. 


2. Two-Way Analysis of Variance (TwoWayANOVA) 

Unlike one-way ANOVA, two-way Analysis of Variance separates the column effects from the 
row effects. Or, in some textbooks, it separates the Treatment effects from the Block effects. 


Example 3: Six air traffic controllers were assigned to test 3 proposals for the modification 
and redesign of their workstations to determine to what extent the 3 alternatives differ in 
terms of ther effects on controller stress. The results for Systems A, B, and C are in columns 
1 , 2, and 3 of the data matrix ate. The result for the six controllers are the rows of the data 
matrix. The Treatments are Systems A, B, and C, the 3 workstation redesigns. The Blocks 
are the six controllers. Since there is expected to be a large variance due to individual con- 
troller differences, two-way analysis of variance is used to control for this variance. 


' 15 

15 

18 ' 

14 

14 

14 

10 

11 

15 

13 

12 

17 

16 

13 

16 

. 13 

13 

13 . 
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( FiT Y 

T 

Ci j 1 h ; ±_ 

Y FHt Y Y FbT V F?tY ) 

IjC- T lr- 1 : i ■■ It - i 1 T .. 4. - 

| 1 001S| 

r E-Way AND'.'A - Block DtsiSn A 

1 g 

ex2_l 

Hnlunn 1 pup! 1 i 


■x2-6 

15 

-List 2: 

ex2_2 


3 

15 

18 

List 3: 

ex2_3 


3 

3 


List 4: 

ex2_4 




List 5: 

ex2_5 




List 6: 

ex2_6 




, CEnter=OK } CESC=CflHCEL^ , 


ex2_6[4]= 

MAIN 

RAD EXACT FUNC fi/ b 



( FlT 

f E-Way AND'.'A - Block D«i3n 


|Tool 

FACTOR: 



ex2_ 

F 

=3. 157895 


15 

P Ualue 

=.057399 


15 

df 

=5.000000 


18 

SS 

=30.000000 



MS 

BLOCK: 

F 

4-P Ualue 

=8.000000 

=5.526318 

=.024181 


ex2 

, CEnter=OK } 

J 


n 

x2_6 

3 

3 

3 


( FlT 

f E-Wa^ AND'.'A - Block D45i3n 

kT 1 


f FlT 


E-Way ANDYA - Block Design 

kT 1 

|Tool 

TBLOCK: 

d 1 


|Tool 

Tdf 

=2.000000 

d 1 

ex2_ 

F =5.526316 

x2_6 


ex2_ 

SS 

=21.000000 

x2_6 

15 

P Oalue =.024181 

3 


15 

MS 

=10.500000 

3 

15 

df =2.000000 

3 


15 



3 

18 

SS =21.000000 

3 


18 

ERROR: 


3 

— 

MS =10.500000 



— 

df 

=10.000000 







SS 

=19.000000 



ERROR: 




MS 

=1.900000 



df =10.000000 








4-SS =19.000000 




s 

=1.378405 


ex2 

^CEnter=0K 5 



ex2 

^CEnter 

O 

II 


MAIN 

RAD EXACT FUNC fi/ b 



MAIN 


RAD EXACT FUNC fi/ b 



I'd like to add another example 111 to illustrate the application - Josef: 

An extra example: Twelve calves are based on their weight at begin of the experiment 
grouped into 4 categories Cl, C2, C3 and C4 with 3 calves in each group. They are fed with 
three different mixtures of feed, FI, F2 and F3. At the end of the experiment the increase of 
weight is measured (see table below). Test the hypothesis that the increase in weigth does 
neither depend on the weight at the begin of the experiment nor on the feed. 



FI 

F2 

F3 

Cl 

7.0 

14.0 

8.5 

C2 

16.0 

15.5 

16.5 

C3 

10.5 

15.0 

9.5 

C4 

13.5 

21.0 

13.5 


#21: animals 


7 

14 

8.5 

16 

15.5 

16.5 

10.5 

15 

9.5 

13.5 

21 

13.5 


#22: TwoWayANOVA(animals) 

Two-Way Analysis of Variance 



Source 

DF 

SS 

MS 

F 

P(F) 


Col umn 

2 

54.125 

27.0625 

5.756277695 

0.0402165846 


Row 

3 

87.72916666 

29.24305555 

6.220088626 

0.02847615714 


Error 

6 

28.20833333 

4.701388888 




. Total 

11 

170.0625 
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#24: TwoWayANOVAC animal s') 


Two-Way Analysis of Variance 



Source 

DF 

ss 

MS 

F 

P(F) 


Col umn 

3 

87.72916666 

29.24305555 

6.220088626 

0.02847615714 


Row 

2 

54.125 

27.0625 

5.756277695 

0.0402165846 


Error 

6 

28.20833333 

4.701388888 




. Total 

11 

170.0625 





The result shows that hypothesis 1 - the starting weight does not influence the increase in 
weight - can be rejected and that hypothesis 2 - the type of feed does not influence the in- 
crease in weight - can be rejected, too. 

3. Two-Factor Analysis of Variance with Interaction (TwoFactorANOVA) 

In some experiments it is desired to draw conclusions about more than one factor or variable. 
In "factorial experiments" all possible combination of factors are involved in the analysis. In a 
two-factor experiment, all combinations include the two factors and their interaction, i.e., fac- 
tor One, factor Two, and One*Two. If there are 'a' levels of factor One, and 'b' levels of fac- 
tor Two, the factorial experiment is considered to have 9 treatment combinations. There 
must be 2 or more samples selected for each treatment combination. 


Example 4: A university wants to determine how three different preparation programs affect 
students' GMAT scores before selecting one to adpot. They include a 3-hour review, a 1-day 
course, and a 10-week course. It is also thought that the college a student is in may affect 
the scores. The university looks at students from the college of business, engineering, and 
arts and sciences. 

There are thus 3 treatments in each of the two factors for a total of 9 treatment cominations. 
The unversity decides to select 2 students for each treatment combination. In experimental 
design terminology, the sample size of 2 for each treatment combination indicates that there 
are 2 replications. 


I'd like to start with the Tl-.procedure, because for applying Don's Derive tool it will be 
necessary to arrange the data in a special way. (I must admit that I had problems to transfer 
the Derive-arranged data into an appropriate form to tackle the problem with the Tl- 
implemented function. I wrote to Don and on the next morning I found his answer in my mail- 
box: 

Josef, 

To use the TI-89 2-Way Analysis of Variance function the data has to be arranged a certain way. 
Factor One becomes the listl, Iist2, Iist3, etc columns in the Stats List/Editor. Factor Two becomes 
the rows of the lists. But the replications, which must be 2 or more, are entered in sequential rows. In 
the following matrix, Factor Two is listed in the rows and Factor One in the columns. 
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Factor Two/Factor One 1 


1 

BOO 

460 

B60 

1 

BSO 

B40 

600 

2 

B40 

B60 

600 

2 

460 

620 

BSO 

3 

4S0 

420 

4S0 

3 

400 

4S0 

410 


I am proud that I nearly had done it. Unfortunately I did not arrange the rows in the correct 
order, I took 1 , 2, 3, 1 , 2, 3 and did not consider the “replications” in the right way. 

I'd like to start with the TI-. procedure, because for applying Don's Derive tool it will be 
necessary to arrange the data in a special way. 


FjT T FHt 

> List Calc 


trl 

tr2 

tr3 

500 

460 

560 

580 

540 

600 

540 

560 

600 

460 

620 

580 

480 

420 

480 

400 

480 

410 





tr3E7]= 


□ istr 


2TT-Test... 

3: 2-SanpZTest... 
4: 2-SanpTTest... 
5: 1-PropZTest... 
6:2-PropZTest... 
7: Chi 2 GOF... 
8:Chi2 2-way... 
9: 2-SanpFTest... 
A: LinRegTTest... 
B: MultRegTests 
C: RHOUfl... 


TYPE DR USE + [EHTER]=DK HMD [ESC]=CAMCEL 


( F1t Y FEt Y Fjt V FHt V F Et Y Fb^ Y F Y 'i 

|Tqq 1 s |P 1 ots |l i si |Ca 1 c |P i sir |Tests 1 1 nts | | 

trl |tr2 |tr3 | -| 1— — ' - 


£-tqqy h ri ■] 1 s i s of Yoriqnct 


Design: 

Luis of Col Factor: 
Lois of Row Factor: 

£Enter=OK 


2 Factor, 

3 + 

B 


Eq Reps-* 


< TESC=CRNCEO 


tr3El]=56G 


TYPE + E E M T E F: ] = □ K liND [E£C]=CAHCEL 


( F1t Y F£t T F]:t Y FHt T F£t T Fb^ T Y 'i 

|Too 1 s |P 1 ots |l i st |Ca 1 c |P i sir |fests 1 1 nts | | 


trl 


50G 

5S0 

540 

460 

4S0 

400 


f £-Wqy ANDYA - £ F qctOK Dtsi3n 

Colum Leuel Lists 


List 1: 

trl 


List 2: 

tr2 


List 3: 

tr3 


, CEnter=0K 5 £ESC=CANCELD ; 


ti*3[l]=56Q 


1 T 


f fit 

Tool 

f £-Wq^ ANDYA - £ FqctOK Dt5i3n 

fi 

COLUMN FACTOR: 

F =1.382872 

P Ualue =.299436 
df =2.000000 

SS =6100.000000 

MS =3050.000000 

ROW FACTOR: 

F =10.269521 

4-P Ualue =.004757 

, £Enter=0K 5 

trl 


500 

580 

540 

460 

480 

400 


tr3 


STAT RAD AUTD FUNC 3/ 3 


r fit 

Tool 


trl 


500 

580 

540 

460 

480 

400 


tr3 


£-Hqy ANDYA - £ FoctQK D*5i3n 


TROW FACTOR: 


F 

P Ualue 
df 
SS 
MS 


=10.269521 

=.004757 

= 2.000000 

=45300.000000 

=22650.000000 


INTERACTION: 

F =1.269521 

4-P Ualue =.350328 

£Enter=OK 


o 


FIt 

Tool 


trl 


500 

580 

540 

460 

480 

400 


tr3 


£-Hqy ANDYA - £ FqctOK D*5i3n 


T INTERACTION: 


F 

P Ualue 
df 
SS 
MS 


=1.269521 

=.350328 

=4.000000 

= 11200.000000 

=2800.000000 


ERROR: 

df =9.000000 

4-SS =19850.000000 

£Enter=OK 




r fit 

Tool 

f £-Wq^ ANDYA - £ FqctOK Dt5i3n ^ 

fl 

Tdf -4.000000 

SS =11200.000000 

MS =2800.000000 

ERROR: 

df =9.000000 

SS =19850.000000 

MS =2205.555556 

s =46.963343 

, CEnter=0K } 

trl 


500 

580 

540 

460 

480 

400 


tr3 


STAT RAD AUTD FUNC 3/ 3 


Interpretation will be given after performing Don's Derive procedure, which is following. 






p 36 
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The results allow us to conclude three things, assuming an a of 5%. First, Factor One, with a 
F-value of 1 .38 and a p-value of .299, does not allow us to reject the Null Hypothesis so we 
conclude that there is no difference in the three test preparation programs on GMAT scores. 
However, Factor Two, with a F-value of 10.269 and a p-value of .00476, allows us to reject 
the Null Hypothesis and conclude that there is a difference in the colleges' ability to prepare 
students for the GMAT. And, F-value and p-value for the interaction effect gives us no rea- 
son to believe that the 3 preparation programs differ in their ability to prepare students from 
different colleges for the GMAT. 


Don wrote: 

Also, for the Two Factor ANOVA I could have programed the function to accept the data 

the same as the TI-89 does. But I chose the way I did because the data does not have to be 
entered in any particular order; the data points could have been entered in any random or- 
der. 


The 3D-plot of a nice Pedal Surface 



pedsurf 




2 2 ' 

> 


2 2 

U + V 



u - 2-v , u + v, 


l 1 

l n 
1 

o 

0 

1 1 

> 

=3 



2 



(see in the next DNL) 
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Surface $ 4i x^y + x' x + y^ z + + 5 z — 0 


********* * * * 



**** ********* * * 
*** ********* ** * * * 
■ * * ************ * 
- *** **************** 
*** • * 


>**•** *** *** *** 

»~* •*•*•*•*•******. 
***************** 
********* ***** 

* * • * *• * •* * • * * 
**•* ***** * * * * 

•aw* 



1 • ***** * *** 
HP— f^J***rT*** ******** * 1 

^♦*** ^***^*7******* ; ] 

«*:*^r****jr/**** **** 

-***£T**** * 1 

^ * * ** * *r IT****** 

*********************** ***** * 



3 3 

VECTOR (Con tourPts_XY(>: -y + z ■>: 


3 3 

+ y ■ z + z + 5 ■ z 


0, z_, -5, 5, -5, 5, 


0.25, 0.25), 


z , -5, 5, 0.25) 


3 3 

ImplicitDotsCx -y + z -x 


3 3 

+ y -z + z +5-z 


0 , [- 5 , - 5 , - 5 ], 


[ 5 , 5 , 5 ], 0 . 2 ) 



Autograph plots 



DPGraph 


Surfer 
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Surface #5: (x 2 + y 2 ) 3 = x 2 y 2 (z 2 + l) 



2 2 3 2 2 2 

VECTOR (Contou rP ts_XY ( (x + y ) = x -y ■ (z +1), z_, -5, 5, -5, 5, 0.25, 0.25), z_, -5, 5, 0.25) 


2 2 3 2 2 2 

VECTOR (Contou rDots_XY((x +y) - x ■y ■ (z +1), 1), 1, -5, 5, 0.5) 



DPGraph Plots 



Surfer Plots (right one with a variation = instead of z 2 +1 simply z 2 ) 
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D. Sjostrand: Incenter and Excenters of a Tetrahedron 
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Incenter and Excenters of a Tetrahedron 

David Sjostrand, Onsala, Sweden 


In this paper we often identify a point P with the vector OP . 

Some definitions 

line(A , B ) is the line passing the points A and B. 

Triangle ( A,B,C) is the triangle with vertices A, B and C. 

Area(A,B,C) is the area of Triangle ( A,B,C) 
d = Area(B , C, A) 
a = Area{C,D,B) 
b = Area(A , C,D) 
c = Area(B,A,D) 

Tetrahedron(A , 5, C, D) is the tetrahedron with vertices ^4, B, C and £>. 
Volume(A , 5, C, D) is the volume of Tetrahedron(A , 5, C, D) . 


Lemma 1 

D divides the line segment BC in the ratio 
(AEC) : (AEB^j counted from C. 

We can immediately verify this lemma using any 
dynamic geometry software. 

Below can see the screens of a Voyage 200 with 
the Cabri-Application: 



a v: Al?i= 1 


$ 

^REE= 1 . 14cm E 
*RCE= 1.27cm* 
\ R: 0.90 

X, BD= 1.70cm 

If \ + DC= 1 . 89cm 

Y^Ox R: 0.90 


. 

T C * 

1 MAIN DEG EXACT FUNC 9 


a I 1 


y ) 

/ 

^REE= 2. 05cm E 
. ^RCE= 2. 90cm E 

\ Rs 0.70 

^X BD= 1 - 4Scm 

\ DC= 2. 10cm 

\ R: 0.70 

In. . . V . " . . ft 


\ THISFDIN-r^^ 



MAIN DEG EXACT FUNC ! 


Proof: 

Area(A,B,D) BD Area(E,B,D) 
Area(A,D,C) DC Area(E,D,C ) 

Area(A,E,B ) BD 

Therefore = . 

Area(A,E,C ) DC 


a I 1 


iRBE= 2. 16cm E 
^RCE= 0. 30cm E 
Rs 2.70 

/ BD=5.70cm 

o/ _^^\\^\_DC=2. 11cm 

NX ^ftJ2.70 

^<5-5 A . , ft 


+ 

MAIN DEG EXACT FUNC j 


Q.E.D. 
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Theorem 1 

The incenter of Tetrahedron(A, B ,C , D) is 
gA ■+■ bB ■+■ cC T dD 
a + b + c + d 

Proof: 

A x e Triangle(B, C, D ) and line(A, A l ) is the line 

having equal distances to the planes 
plane(A, B, C ) , plane(A, B, D ) and 

plane(A,C,D ) . 


A 



We define points B t e Triangle(A, C, D ) , C x e Triangle(B, A, D ) and D x e Triangle(B, C, A) and 
lines line(B, B x ) , line(C, C \ ) and line(D, /), ) in a corresponding manner. The incenter of 
Tetrahedron(A, B, C, D ) is the intersection point of these four l in es. 

Let e be the common distance between A x and plane(A, C, D ) and between A x and 
plane(A,B,D ) . 

be ce 

Now Volume(A, A X ,C,D) = — and Volume(A, A X ,B,D) = — . 


Let h a be the altidude of Tetrahedron(A, B, C, D ) from A to plane(B, C, D ) . We receive 

. , . , „ ^ h area(A,,C,D) . , „ „ „ _ h area(A,,B,D) 

Volume(A, A x ,C,D) = — and Volume(A, A x ,B,D) = — ! . 


This means that 


be h a area(A x ,C,D ) 

b _ 3 _Volume(A,A x ,C,D) _ 3 _ area(A x ,C,D) 

c ce Volume(A,A x ,B,D ) h^area^A^B^^ area(A x ,B,D ) 

3 3 


By Lemma 1 we now conclude that A',, divides the line segment SC in the ratio b:c counted from C to 

bB + cC 


B. This means that X n = ■ 


b + c 


By symmetry we have that X c = 
(b + c\X„ +dD 

The point 
The point 
Therefore 


bB + dD 


b + d 

(b + c)X D +dD bB + cC + dD . , v ^ 

= is on LineiX n ,D) 

b+c+d b+c+d D 

(b + d)X c +cC _ bB + dD + cC _ bB + cC + dD . 


{b + d^ + c 
bB + cC + dD 


is on line{X c , D ) . 


4 = 


b + c + d 
bB + cC + dD 


b + d + c b + c + d 

is the intersection point between line(X d ,D ) and thus 


b + c + d 




D-N-L#69 D. Sjostrand: Incenter and Excenters of a Tetrahedron p 41 


By symmetry B x = 


aA + cC + dD 
a + c + d 


(b + c + d \ A x + + bB + cC + dD . 1 . / . . x 

As A — = is on line line(A , 

(& + c + £/j + $ a+b+c+d 


( + c + d ^ ^ + bB aA + bB + cC7 + dD . 


((2 + c + ^/^ + Z? cl b c d 


is on line line(B,B x ) , we see that 


£l4 + bB + cC + dD 
cl -\- b -\- c -\- d 


is the intersection point between line(A , and line(B , 5,) . 


Thus the incenter of Tetrahedron(A,B,C,D ) is 
Q.E.D. 


+ ZrZ? + cC + dD 
cl -\- b -\- c -\- d 


Theorem 2 

The excenters of Tetrahedron(A,B,C,D ) are given by 

-aA + bB + cC + dD aA-bB + cC + dD aA + bB -cC + dD aA + bB + cC-dD 

, , and . 

-a+b+c+d a-b+c+d a+b-c+d a+b+c-d 

We will first prove that the excenter of Tetrahedron(A,B,C,D ) lying opposite to the plane defined 


by A, D and C counted from 5 is given by 


aA-bB + cC + dD 
a-b + c + d 


Proof: Everything we did in the proof of Theorem 1 can be repeated word by word with the only ex- 
ceptions that we have to replace b by -b in the expressions for X c and X D and in all the subsequent 
expressions. It is straightforward to justify this. 


If the line having equal distances to the 
planes plane(A,B,C ) , plane(A,B,D ) and 
plane(A , C, D) is parallel to 
plane(B,C,D ) we argue as follows: 


aA-bB + cC + dD 


and excenter of Tetrahedron(A,B , C,D) lying opposite to the plane 


a-b+c+d 

defined by A, D and C counted from 5 both are continuous functions of the vertices of the tetrahedron 


it follows that the expression, 


aA-bB + cC + dD 
a-b + c + d 


for the excenter, is valid even in this case. 


By symmetry the three other excenters of Tetrahedron(A,B , C,D) are given by 

— aA + bB + cC + dD aA + bB — cC + dD aA + bB + cC — dD 

, and . 

-a+b+c+d a+b-c+d a+b+c-d 


Q.E.D. 





p 42 


D. Sjostrand: Incenter and Excenters of a Tetrahedron 


D-N-L#69 


This is the /)/T<7r/>application of the results from above together with a graphic representation of the 
tetrahedron together with the spheres: 

Inspheres and Exspheres of a Tetrahedron 

David Sjostrand 

#1: CaseMode := Sensitive 

#2: InputMode := Word 

Area of a triangle with vertices A, B and C: 

1 

#3: area(A, B, C) := | (B - A) x (C - A) | 

2 

Area of the faces of a tetrahedron with vertices A, B, C and D: 

#4: [a := area(B, C, D), b := area(A, C, D), c := area(B, A, D), d := area(A, B, C)] 

Volume of a tetrahedron with vertices A, B, C and D. 

Radius of the inscribed circle of this tetraedron. 

l 3 • v 

V := | (D - A) • (B - A) x (C - A) | , radius(a, b, c, d) := 

6 a + b + c + d. 

Sphere with radius r and center M: 

Sphere(r, s, t, M) := |"M + r-COS(s) -SIN(t) , M + r • SIN(s) • SIN(t) , M + r-COS(t)"| 

#6: L 1 2 3 J 

Incenter of a tetrahedron with vertices A, B, C and D: 

l 

#7: incenter(a, b, c, d) := (a-A + b-B + c-C + d-D) 

a + b + c + d 

#8: pi ane(X, Y, Z) DET([X - Y, X - Z, [x, y, z] - X]) = 0 

#9: [plane(A, B, C) , plane(A, B, D) , plane(A, D, C) , plane(D, B, C)] 

The insphere: 

#10: Sphere(radi us(a, b, c, d) , s, t, incenter(a, b, c, d)) 

The four exspheres: 


#11: 

Sphere(radi us(-a, b, c, 

d), s, 

t, 

i ncenter(-a 

, b, 

C, d)) 



#12: 

Sphere(radi us(a, -b, c, 

d), s, 

t, 

i ncenter(a, 

-b, 

c, d)) 



#13: 

Sphere(radi us(a, b, -c, 

d), s, 

t, 

i ncenter(a, 

b, 

-c, d)) 



#14: 

Sphere(radi us(a, b, c, 

-d), s, 

t, 

i ncenter(a, 

b, 

c, -d)) 



#15: 

[A := A, B :: B, C := C, D 

:= D] 







#16: 

i — i 
> 

l \ 

i — i 

o 

o 

0 

TO 

1 f 

i — i 

I— 1 

I 1 

o 

o 

C : 

1 1 

i — i 

o 

i— 1 

0 

1 1 

D 

i — i 

O 

O 

1]] 


#17: 

[A [-3, -4, -5], B := 

[1, 2, 

0], 

C [0, 1, 

3], 

D := [0, 

-2, 

1]] 

#18: 

1 1 

> 

*1 \ 

l 1 

o 

0 

1 

1 1 

00 
*1 \ 

1 — 1 

, 4, 3] 

, c 

:= [-4, 5, 3 

], D 

:= [5, - 

-5, 

6]] 
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Here are the four planes together with the spheres. The insphere is invisible. 



The next plots show the tetrahedron with its edges together with the four spheres. In 
the second plot I changed the colour and didn't show the mesh lines. 


#19: [A, B, 

, c, 

D 

, B, 

D, 

A, 

C] 






" A 

B " 


' A 

D " 


' B 

D " 


' A 

c r 

#20: 

A 

C 

■ 

A 

B 

■ 

B 

C 

l 

A 

D 




Expression #19 plots the wire grid of the tetrahedron. 

Expression #20 plots the sides of the tetrahedron. 

David asked, if I could find a proof Lemma 1 for an arbitrary position of point E (even outside of the 
triangle). I answered that in my opinion his proof is also valid for these positions but I felt inspired to do 
another proof which is also valid for all positions of point E. 
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Appendix: Another proof of Lemma 1 : 

#1: [InputMode := Word, CaseMode := Sensitive] 

e2 - a2 

#2 : y - a2 (x - al) 

el - al 

e2 - a2 

#3: 0 - a2 = (dl - al) 

el - al 

' e2 - a2 

#4: SOLVE 0 - a2 (dl - al) , dl 

«» el - al 

al*e2 - a2*el 

#5: dl = 

e2 - a2 

al-e2 - a2-el 

#6: BD := 

e2 - a2 

#7: DC := BD - cl 

'XX 1 
1 2 

1 Y Y 1 

#8: area(X, Y, Z) := DET 1 2 

2 

Z Z 1 
. 1 2 

#9: [A := [al, a2], B := [0, 0], C := [cl, 0], E := [el, e2]] 

area(A, E, B) 

area(A, E, C) 

#10: = 1 

BD 

DC 

Q. E. D. Now in detail: 

al-e2 - a2-el 

#11: area(A, E, B) = 

2 

al-e2 + a2*(cl - el) - cl*e2 

#12: area(A, E, C) = 

2 

area(A, E, B) al-e2 - a2-el 

#13: = 

area(A, E, C) al*e2 + a2*(cl - el) - cl*e2 

BD al*e2 - a2*el 

#14: = 

DC al-e2 + a2-(cl - el) - cl-e2 


David's comment: Your dfw-file is interesting. As I understand it, the four subexpressions of 
area(A, E, B)/area(A, E, C)/(BD/DC) can be negative or positive, which seems to be nice. 






